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Abstract. We investigate the mesoscopic disorder induced rms conductance variance 
SG in a few layer graphene nanoribbon (FGNR) contacted by two superconducting 
(S) Ti/Al contacts. By sweeping the back-gate voltage, we observe pronounced 
conductance fluctuations superimposed on a linear background of the two terminal 
conductance G. The linear gate-voltage induced response can be modeled by a set 
of interlayer and intralayer capacitances. SG depends on temperature T and source- 
drain voltage V s d- SG increases with decreasing T and \V s d\- When lowering |V^d|, a 
pronounced cross-over at a voltage corresponding to the superconducting energy gap 
A is observed. For \V s <i \ $ A the fluctuations are markedly enhanced. Expressed in the 
conductance variance Ggs of one graphene-superconductor (G-S) interface, values of 
0.58 e 2 /h are obtained at the base temperature of 230 mK. The conductance variance 
in the sub-gap region are larger by up to a factor of 1.4 — 1.8 compared to the normal 
state. The observed strong enhancement is due to phase coherent charge transfer 
caused by Andreev reflection at the nanoribbon-superconductor interface. 



PACS numbers: 72.80.Vp, 73.23.-b, 73.40.-c, 74.45. +c 
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1. Introduction 

Graphene nanostructures [U [2] provide the unique opportunity to study fundamentally 
new quantum coherent phenomena and, on the other hand, possess an immense potential 
for applications [3]. The new quantum phenomena in graphene originate from the 
linear band structure of the relativistic-like quasiparticles and their chirality. The 
demonstration of a tunable supercurrent posed fundamentally new questions regarding 
the sources of apparent dephasing in this material since the amplitude of the critical 
current is lower than expected [U [5] . Previous studies in disordered graphene devices 
have shown [6j[7J[8] that interference effects such as the universal conductance fluctuation 
(UCF) and weak localisation are dominant corrections to the conductance at low 
temperatures. First studies of UCF and weak localisation were done in the 80 's in 
metallic systems and 2D semiconducting heterostructures [9]. Conductance fluctuations 
arising in these systems are universal, independent of sample size and degree of disorder, 
reaching values of the order e 2 /h [9]. The interest in studying conductance fluctuations 
in graphene stems from predictions that they cease to be universal in the coherent 
state of disordered graphene [10]. In graphene samples a range of sources of disorder 
has been identified: the close proximity with the substrate [TT], interaction with the 
leads [El [13], rippling of the graphene layer [13], and unintentional doping. The disorder 
reduces the mobility of the quasiparticles. In the strong disorder regime, the variance of 
the conductance SG coincides with the predicted value for disordered metals, whereas 
in the weak disorder regime, 5G is larger than the universal value due to the absence 
of back-scattering, characteristic of the honeycomb lattice of graphene [10]. Additional 
information on UCF in the system can be gained by attaching superconducting contacts 
to the disordered region |15j . At the interface between a normal metal N and a super- 
conducting electrode S phase sensitive Andreev reflections occur for energies lower than 
the superconducting energy gap A. At the N-S interface, an electron coming from N 
couples with its time-reversed counterpart to form a Cooper pair which can enter into 
S. The phase coherent Andreev states at the N-S interface can be destroyed by applying 
a magnetic field which breaks time reversal symmetry, by increasing the temperature, 
or by applying a large source-drain voltage. 

In this work we study the two terminal conductance of a few layer graphene 
nanoribbon (FLGNR) contacted by Ti/Al leads as a function of back-gate voltage V g 
and source drain voltage V s d- We find that the conductance variance increases with 
lowering temperature reaching an amplitude of the order of e 2 /h at zero bias and 
base temperature of 230 mK. By applying a source drain voltage larger than A/e 
the conductance variance decrease by up to a factor of 1.8. A characteristic cross- 
over at an energy corresponding to A confirms that the observed enhancement is due to 
Andreev reflection at the graphene-S (G-S) interface. This finding complements existing 
work [3 [8] by focusing at the bias dependence of conductance variance. 
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2. Sample fabrication and characterisation 

Graphene flakes are prepared by mechanical exfoliation of natural graphite (NGS GmbH, 
Leinburg, Germany) using a surface protection tape (SPV 224P from Nitto Denko), 
followed by the transfer of flakes onto a piece of highly p-doped Si wafer with a top 
thermal oxide layer of thickness t = 304 nm. The high-doping of the substrate ensure 
the possibility to gate the flakes by applying a back-gate voltage V g to the substrate. 
After flake transfer the samples are rinsed in solvents to remove the glue residue from the 
flakes and substrate. Suitable flakes are selected and localized by an optical microscope 
with respect to a grid of markers. The devices were patterned with e-beam lithography 
and subsequently metallized in an UHV e-beam evaporation system at a pressure of 
10 _7 mbar, followed by lift-off in acetone. We focused on single and few-layer graphene 
devices contacted with superconducting aluminium (Al). In particular, we studied the 
sample shown in the inset of figure 1. This narrow multi-layer flake was contacted with 
a Ti/Al/Ti tri-layer (5/40/20 nm). The purpose of the bottom Ti layer is to ensure 
high contact transparency between Al and graphene, whereas the top Ti layer caps the 
Al underneath. From the SEM image, we estimated the width W of the flake to be 
between 150 nm and 200 nm and the edge-to-edge distance L between the leads to be 
~ 225 nm. 

The back-gate dependence of the two terminal linear differential conductance 
G = dl/dV taken at 1.4 K and 5.5 K (shifted by 0.2 • 2e 2 /h for clarity) in the range 
of V g G [— 10 V, 50V] is shown in figure 1. Both traces are taken at temperatures where 
the contact leads are in the normal state. They show a minimum of the conductance 
Gmin ~ 3.2e 2 /h at 20 V, marking the position of the charge neutrality point (CNP). 
A small back-gate shift between the two curves of ~ 2 V is observed after cooling the 
sample from 5.5 K to 1.4 K. At 5.5 K the conductance linearly depends on V g away from 
the CNP and shows small aperiodic, but reproducible fluctuations. The amplitude of 
the fluctuations becomes larger after lowering the temperature to 1.4 K. 

The magnitude of the conductance slope AG/AV g (marked with dashed lines in 
figure 1) can be used to determine the number of layers as demonstrated in the work of 
Zhang et al. [16]. Because of the strong electrostatic screening in the vertical direction 
of the FLGNR stack, the outermost carbon layer is affected most by the back-gate 
voltage. In the inner layers, in contrast, the electric gate- field is strongly suppressed so 
that these layers add a nearly gate-independent conductance value proportional to the 
number of layers to the total conductance G. If one assumes that the electron diffusion 
constant is the same in all layers, an estimate for the number of layers Nl in the FLGNR 
can be obtained from the change in conductance AG/AV g normalized to the minimum 
conductance G min [TH E] ■ 

To be more explicit, we use the two-dimensional model described in figure 2a. The 
back-gate is spaced by a distance t from an infinite set of graphene layers j = 1,2,3,.... 
The interlayer thickness between the graphene layers is taken to be d. Ej, Qj, and 
4>j denote the electrical field, the areal density of the excess charge measured from the 
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Figure 1. Two terminal linear conductance G taken at 5.5 K and 1.4 K as a function 
of the back-gate voltage V g . The 5.5 K curve is shifted by 0.2 • 2e 2 /h to clearly show 
the increase of the conductance fluctuations upon lowering the temperature. The inset 
shows a scanning electron micrograph of the sample with Ti/Al/Ti electrodes across 
a few-layer graphene nanoribbon. 

CNP, and the electrostatic potential in the different graphene layers (j > 1) or back- 
gate (j=0), respectively. For simplicity we set the dielectric constant to the vacuum 
permittivity eo- This assumption can be relaxed afterwards by replacing the relevant 
parameter by the correct back-gate capacitance C g . 

The difference in the electric fields is determined by the excess charge according to 

Q; = 60^ (1) 

We assume that the graphene stack has a back contact (in our case realized by the 
source and drain contacts), which are set to zero, whereas the back-gate is biased to an 
electrochemical potential V g . Because electrons can be exchanged between the different 
graphene layers by tunnelling, the electrochemical potentials of all layers are equal and 
zero in the thermodynamic limit. As a consequence the chemical potential fij of layer j 
is given by the negative of the electrostatic potential <pj, i.e. \ij = —ecfij with e the unit 
of charge. On the other hand, the chemical potential \Xj is determined by the excess 
charge density according to Qj/eNj, where the new symbol Nj denotes the areal density 
of states in graphene layer j. Taken together, we arrive at 

0. = -Q j/e "N r (2) 
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Figure 2. (a) Schematics of the model to calculate the electric-field screening in stack 
of few layer graphene which is gated by an electrostatic back-gate voltage V g . Ej and 4>j 
denote the electric field and the electrostatic potential in layer j. (b) The circuit can be 
reduced to a series connection of two capacitances, where Cq is an effective quantum 
capacitance and C g the gate capacitance, (b) The full circuit diagram where C g , 
Ci, and C q denote the gate, graphite interlayer, and graphene intralayer capacitance, 
respectively. The latter capacitance is sometimes also termed quantum capacitance. 



Adding for the charge Qj the values given by (CQ), yields 
_ e (£j-i - Ej) 



em, 



(3) 



A self-consistency relation can now be formulated by noting that the difference of the 
electrostatic potentials determines the electrical fields, i.e. <pj — <j)j-i = Ejd. This leads 
to 

where the density-of-state has been assumed to be constant and given be N. This 
equation has exponentially decaying solutions of the form 

Ej = E e~ d ^\ (5) 

where A denotes the interlayer screening length. Placing (jSj) into (jl]) yields the following 
condition: cosh(d/X) — 1 = e 2 Nd/2eo- This equation determines the screening length 
A. In order to quantify A, the areal (2d) density-of-states N is required. We estimate 
N from the known 3d value N of graphite to iV = N d. With N = 5.2 ■ 10 20 cm~ 3 /eV 
and the interlayer distance <i = 0.34nm [18], one surprisingly obtains d/X = 1, i.e. 
A = 0.34 nm. Eventually, we can place the solutions Ej in (jSJ) back into the equations 
(T5]) and ([3]) to obtain the charge in layer j: 

Qj = -Qoe-^'- 1 ) {l - e- d ' x ) . (6) 

Obviously, if screening is strong (A << d), \Q\\ ~ \Qo\, whereas |Qj>i| << \Qo\- In the 
opposite case of weak screening (A >> d), the charge decays slowly and is given for not 
too large j by \Qj\ ~ \Qo\d/ A. In the present situation with A = d, 63 % of the charge is 
in the first layer, 23 % in the second, and the third layer already carries less than 10 %. 
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Hence, this model of an infinite stack of layers should work well for few layer graphene 
stacks when the number of layers Nl ~ 3. The parameter Q in ([6]) can be related to 
the applied gate voltage V g using the relation (fio — (pi = V g — <fii = E t. We then arrive 
at a relation between the charge Q on the gate and the gate voltage V g : 

Qo = i/Cq + iav (7) 

where C g is the geometrical gate capacitance and Cq is an effective chemical capacitance 
(sometimes also termed quantum capacitance), both taken per unit area. Cq is given 
by 

e 2 N 

C Q = i- e -d/x- ( 8 ) 

Equation [7| shows that this relation represents a series connection of C g with an effective 
quantum capacitance Cq as depicted in figure 2b. The latter can also be derived from 
the circuit shown in figure 2c consisting of an infinite series of interlayer capacitances 
Ci and intralayer quantum capacitances C„ = e 2 N [19]. We next estimate the two 



capacitances C g and C q . Taking e = 3.9 for the relative permittivity of Si0 2 yields 
C g = 1.1 • 10 _4 F/m 2 . With the aerial density of states N, which we estimated before 
from known graphite values, we obtain C q = 2.8 • 10 _2 F/m 2 . Hence, C q ^> C g by more 
than two orders of magnitude. Since the smaller capacitance counts in a series connection 
as the one shown in figure 2b, the relation between gate charge Qo and gate-voltage V g 
is to a very good accuracy given by the normal one Qo — C g V g . 

Having analysed the screening problem, we can calculate the sheet conductivity of 
the whole stack. We assume to be in the diffusive limit and use the Einstein equation, 
which relates the conductivity aj to the density-of-state in layer j and the diffusion 
constant D. In a Fermi gas the diffusion coefficient is given by D = vpr/2, where 
vp is the Fermi velocity, in graphene equal to 10 6 m/s, and r the scattering time. For 
simplicity we assume that D is constant. This is an approximation, as it is known 
that there are surface effects and it is plausible that adsorbates must have the biggest 
effect on the mobility of the first layer. With this assumption we can write for the total 
conductance G of the FLGNR 
W 00 

G^y-e^Ety, (9) 

where W and L are the width and length of FLGNR, respectively. The gate-dependence 
of the conductance is due to the energy dependence in the density-of- states Nj(E) 
which we have to add now here. In an ideal single graphene layer the energy dispersion 
relation is given by E = hvp\k\, where k denotes the wavevector in two dimension. This 
dispersion results in a density-of-state N given by [H] 



n ^ = ? < io) 
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This density-of-state goes to zero for E — ► 0. We have however assumed a finite density- 
of state N d at the charge neutrality point (CNP, i.e. E — 0), which is caused by the 
interlayer overlap. To interpolate between the two regimes, we write: 

N{E) = N d + -tJ^L = AT rf(l + p\E\). (11) 

The slope (3 is not a free parameter, but determined by the above equation. We obtain 
(3 = 8.3 /eV. Since the energy E in (fTPl denotes the chemical potential, we can replace 
it with the electrostatic one. This leads together with Qj to: 

W °° 

G = e ^ DNod ^2(l + /3|e&|). (12) 

Adding the explicit expressions for 4>j(V g ), the result is 

^^(^(..^-^j. (!3) 

By dividing with the minimum conductance G min at V g = 0, the dependence on the 
diffusion coefficient drops out. In the practical limit of C q ^> C g the gate voltage change 
of G/G min is given by the final simple result: 

d ( G \ _ f3e C g 



, , , • (14) 

dV g \G min J N L C q 

We can apply this result to figure 2. Equation [TJ] predicts a relative change of 3.3 %/V 
for Nl = 1. We measure a change of 1.2 %/V on the right and 0.32 %/V on the left. 
These two slopes correspond to Nl = 3 on the right and Nl = 10 on the left. The model 
clearly shows that we are dealing with a number of layers, but that we are in the regime 
of few layers with Nl < 10. Due to the different slopes for voltages smaller or larger than 
the CNP, a more accurate estimate for Nl cannot be given. However, we stress that 
physically we deal with one FLGNR with one given Nl number. The reduced change 
on the gate-voltage to the left of the CNP suggests that the carrier density is markedly 
increased on the hole side leading to an enhanced quantum capacitance. As the CNP is 
strongly on the positive side, the FLGNR is substantially hole doped. This hole doping 
may be caused by states that are induced by the source and drain contacts. 

The dependence of conductance G(V g ) on the gate voltage V g can also be used to 
determine the mobility fi. Taking the derivative of equation ( ITBl versus V g and using 
the relation C q ^> C g yields: 

W g = T De ^- (15) 

With the mobility fi one may also write the conductance as 
W 

G = —nQ, (16) 

where Q is the effective carrier density. Since C q ^> C g , the gate- voltage induced carrier 
density is given by Q = C g V g = Qq. This leads to 
dG W „ 

W a = T" c »' (17) 
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By comparing this equation with equation (fT5|) . we deduce the relation \i = Dej3. Taking 
the experimental value for AG/AV g ~ 1.2e 2 /h per 30 V, equation ffl7|) provides the 
value /i = 140cm 2 /Vs for the electron mobility. D follows then to be ~ 17cm 2 /s, from 
which we estimate by virtue of D = vpl e /2 the scattering mean-free path to l e ~ o . o urn. 
With l e < L, W, we conclude that the device is diffusive as anticipated in the beginning. 

3. Conductance measurements 

Upon cooling the sample below the critical temperature of bulk Al, = 1.2 K, an 
energy gap A bulk pa 200 /ieV opens in the density of states of Al. At the base temperature 
of 230 mK the transport properties of the device are deduced by measuring the two- 
terminal differential conductance G as a function of source-drain V s d and back-gate 
voltage V g , represented by the grey-scale plot in figure 3a. The differential conductance 
G was obtained by superimposing an AC excitation voltage of 20 /iV onto the DC part 
of the source-drain voltage and measuring the pre-amplified AC current with a lock-in 
amplifier. For each back gate voltage in the range from —10 V to 20 V V s d is swept from 
-0.8 mV to 0.8 mV. 

In the grey-scale plot shown in figure 3a we observe a reproducible pattern of 
conductance fluctuations which is modulated by both V g and V s d- The amplitude of 
conductance fluctuations clearly diminishes as V s d is increased. In contrast, as a function 
of back-gate voltage the fluctuations have a homogeneous amplitude, similar to what 
is seen in figure 1. As compared to the data in figure 1, the amplitude is larger in 
figure 3a due to the lower temperature. This dependence of conductance fluctuations 
points to so-called 'universal conductance fluctuations' (UCF), which are caused by 
quantum interference effects [9]. 

In addition to UCF, the presence of superconductivity is evident from the 
pronounced dark (low) conductance band around zero bias in between the 2A-lines, 
which are marked in figure 3a by white dashed lines. This band disappears at 
temperatures above the critical temperature T c , as well as in a magnetic field larger 
than the critical field. The band of reduced conductance shows, that the conductance 
is on average lower at low bias than it is at high bias (normal state). We illustrate this 
by overlaying the average conductance G(V s d) versus V s d in the middle of the grey-scale 
plot. This, however, does not mean that all individual G{V s d) curves show dips at zero 
source-drain voltage. Though the majority of G(V s d) curves must display dips, such 
as the example overlayed on the left side of the grey-scale plot, there is a significant 
number of curves that display peaks around zero bias instead. An example is shown by 
the right overlayed curve (see also figure 4). 

To quantify the energy scale at which the superconductivity related effects become 
important, we extract the critical temperature T c in the following. The conductance 
G(T) as a function of temperature T was measured during the cool down at V g — 0. 
This is shown in figure 3b. For this gate- voltage value, G(T) decreases with temperature, 
reaching a minimum at about 0.6 K. Below this temperature the conductance sharply 
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Figure 3. a) Grey-scale plot of the two-terminal conductance G = dl/dV as a 
function of the back-gate voltage V g and the source drain voltage V s d at T = 230mK. 
The dashed white lines indicate the position of twice the superconducting gap 2A/e. 
An example of a conductance spectra with dip around zero bias is shown on the left, 
whereas one with a peak is shown on the right with the average curve placed in the 
middle of the plot, b) Temperature dependence of the linear conductance taken during 
the sample cool-down at V g = 0. (c) The conductance averaged over all back-gate 
curves G{V s d) as a function of V s d- (d) In our system the interface resistances between 
the graphene ad superconductor dominate over the internal resistance of the FLGNR. 
This is illustrated in the resistor network. 



increased reaching a value of 1.72- 2e 2 / h at the base temperature. We take the cross-over 
temperature at 0.6 K as the transition temperature T c to the superconducting state. We 
note that the bulk value for Al is T c = 1.2 K. The substantial reduction of T c of 0.6K in 
our superconductor-graphene-superconducutor (S-G-S) device as compared to the bulk 
value may be due to the inverse proximity effect from the two normal metal Ti layers 
with thicknesses of 5nm and 20 nm surrounding the Al [20J. 

Figure 3c shows the average conductance G(V s d) as a function of V 8 d- This curve is 
obtained as the mean of all individual conductance curves measured at a fixed V g value, 
after a small linear background was subtracted from the grey-scale plot in figure 3a. In 
order to fit the data to the BTK model [21] , we need to know how the voltage drops in 
our S-G-S device. To estimate this, we evaluate the resistor model shown in figure 3d. 
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The total resistance is divided into two interface-graphene resistors Rgs, which on both 
sides are taken to be equal, and the intrinsic resistance Rq of the FLGNR connected 
in series. We compare the measured minimum conductance G m i n with the expected 
minimum conductivity a min of a single sheet of graphene. To do so, we will use the 
number of graphene layers in our FLGNR which we estimated in the previous section 
to be between 3 and 10. 

Unlike the electrons in a two-dimensional electron gas with quadratic dispersion, the 
Dirac-electrons do not localize, even not for the smallest possible carrier density. Instead, 
there is a minimum conductivity a m i n which is reached at the CNP. In experiments with 
graphene layers of intermediate disorder, a m i n has been found to be ~ Ae 2 /h [22J. In 
the ballistic regime a m i n (4/ir) e 2 /h has been predicted [23]. This value, which is reduced 
by a factor of it, has not yet been demonstrated experimentally, even not in the highest 
quality suspended graphene devices [21]. 

The estimate of the elastic mean-free path in the previous chapter has shown that 
we are in the diffusive limit. In our experiment G m i n ~ 3.2e 2 /h, which relates to 
< le 2 /h per layer. Since this value is appreciably smaller (by ~ a factor of 4) than 
a m i n in the diffusive limit [22] , the total resistance must be dominated by the interface 
resistances, i.e. IRqs > Rg- Taking this into account, we simplify our network by 
dropping R G altogether. We therefore assume that the applied source-drain voltage V s d 
drops symmetrically over the source and drain contact. The voltage drop Vgs across 
one interface is then half of the voltage V S( i applied across the junction and, hence, 
Ggs = 2(7. 

We are now in the position to model the measured average conductance curve 
G(V s d) in figure 3c using the BTK model [25l l21~] . This fit, which is shown in 
figure 3c by the dashed curve, provides as one parameter the superconducting energy 
gap A av = 60 /ieV. Other fitting parameters deduced from the model are the effective 
barrier strength Z « 0.6, and an inelastic broadening parameter r ~ 10 fieV. The 
broadening parameter is introduced as an imaginary energy term in the BCS density 
of states, p(E,T) — (E — iT)/[(E — iT) — A 2 ] 1//2 with E being the energy of quasi- 
particles measured from the Fermi energy [25]. The superconducting energy gap A av 
can be compared to the previously deduced transition temperature T c = 0.6 K. Taking 
the BCS relation A = 1.76/jbT c [20], one obtains A « 90 fieV. Good agreement between 
the two numbers is found. The consistent numbers support our simplified approximation 
that most of the voltage drops over the superconductor-graphene contacts. 

Figure 4 provides examples of individual conductance curves G(Vgs) taken for 
different values of V g . Note, Vqs is the voltage over one graphene-superconducutor 
junction and taken to be half of V s d- As we have emphasised before, G(Vas) is suppressed 
on average in the sub-gap region for small source-drain voltages V s d because the majority 
of curves show dips around zero V s d- Examples of such curves are shown in figure 4a for 
back-gate values V g = —7.7 V, —6.2 V, 16.5 V and 9.9 V. However, a substantial number 
of curves show peaks. Examples are given in figure 4b taken at V g — —5.8 V, 8.3 V, 
8.5 V and 11 V. This shows that Andreev reflection can dominate sub-gap transport 
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Figure 4. Typical differential conductance curves G(Vqs) for different constant back- 
gate voltages V g displayed as a function of the voltage Vqs that drops over a single 
graphene-superconductor junction. Vqs — V s d/2- (a) In the bias window \Vgs\ < A/e 
a dip is observed in G(Vgs) for V g = —7.7, —6.2, 16.5, and 9.9 V, whereas (b) a peak 
appears for V g = —5.8, 8.3, 8. 5, and 11V. All curves are normalized to the large bias 
conductance values G N which are 1.76 (-7.7), 1.85 (-6.2), 2.44 (16.5), 2.38 (9.9), 1.92 
(-5.8), 2.23 (8.3), 2.25 (8.5), 2.48 (11) in units of 2e 2 /h, where the number in bracket 
denotes V g . 

at the G-S interface and thereby enhance G. After fitting individual curves to the 
BTK model [251 [21] (dashed lines in (a) and (b)), we extract a superconducting gap of 
A/e ~ 65 /xV, which agrees very well with A av = 60 /xeV deduced for the average curve 
before. The inelastic broadening parameter T ranges between 10 — 20 /xV and for the 
barrier strength Z we obtain « 1 for curves with dips and Z ~ 0.3 for traces with peaks. 
This agrees with the notion that the conductance of a superconductor-normal interface 
can be enhanced in the superconducting state, provided the interface transmission is 
sufficiently high, and hence, the barrier strength Z is sufficiently small. 

We now turn our attention to the bias-dependence of the conductance fluctuations. 
From the conductance data in figure 3 obtained at 230 mK we determine the rms 
conductance variance for a fixed source-drain voltage V s d by a statistical average over all 
G(V g ) values in the gate-voltage range. We assume that the large gate-voltage window 
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provides a sufficiently large ensemble of different disorder configurations. In order to 
determine the root-mean square variance 8G a, linear background has been subtracted 
from G{Vg) leaving us with values AG(V g ) with a zero mean. The rms variance 5G of 
the whole two-terminal device is then given by (AG 2 ) 1//2 . For comparison with theory it 
can be convenient to relate the conductance variance to one graphene-superconductor 
(G-S) interface, which we denote by SGqs- If we assume that the two G-S contacts 
fluctuate independently 6Gqs — V2(5G). This is plotted in figure 5 as a function of 
voltage Vqs — Vsd/2 dropping over one contact. 




0.00 0.05 0.10 0.15 0.20 0.25 0.30 



V GS (mV) 

Figure 5. The conductance variance of a single graphene-superconductor (G- 
S) interface SGqs (root mean square of the deviation of G(V s d) relative to linear 
background times y/2) as a function of the voltage drop across the G-S interface. 
8G(Vqs) displays a pronounced increase below a cross-over voltage that agrees with 
the superconducting gap value A. The inset displays G(V g ) measured at different 
temperatures 5.5 K, 1.4 K and 0.23 K, after subtracting a linear background for 
comparison. The 5.5K curve is shifted by 0.2-2e 2 //i and the 0.23K one by —1.2-2e 2 /h 
for clarity. 

In figure 5 we note two distinct regions with a clear cross-over at the energy 
eVos ~ 90/ieV (arrow), which coincides with the previously estimated superconducting 
energy gap A of the Ti/Al/Ti tri-layer. In both regions SGqs decreases with the applied 
bias. However, the decreasing rates of SGqs are different for the two regions. In the 
high bias region, Vqs > A/e, 5Ggs decreases four times slower than in the region below 
the superconducting energy gap A. In the latter region, SGqs starts saturating as the 
voltage drop across the junction Vgs approaches fc^T/e ~ 20 /iV and takes a finite value 
of SGos,230mK = 0.58e 2 //i at zero bias. In addition, from the inset where conductance 
traces at 5.5 K, 1.4 K and 230 mK are shown, we find for the conductance variance at 
1.4 K the value SGgn,iak = 0.19e 2 //i and at 5.5 K SGgn,5.5K = 0.06 e 2 /h, where the 
index N denotes that Al is in the normal state. 
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It is clear that the observed increase in SGgs is due to the presence of the super- 
conductor and Andreev reflections at the G-S interface. At zero temperature and zero 
magnetic field theory predicts an increase of 2.07 for the conductance variance of a single 
superconducting-normal metal junction in the superconducting relative to the normal 
state [15] . A simple estimate provides from our data an enhancement of ~ 1.4, which is 
smaller than the maximum that can be expected. This estimate is obtained by linearly 
extrapolating SG from the high bias down to zero bias (lower dashed line in figure 5). 
This extrapolation yields SG = 0.3e 2 /h for the variance SG of the whole two-terminal 
sample. 

The reduced magnitude of the conductance variance ratio SGgs/SGgn is not 
unexpected, taking the energy of the lowest temperature and the finite AC-bias 
modulation into account, which both lead to saturation of SGgs at the lowest bias 
voltage in figure 5. The reduction of the saturation value 5G from the expected 
'universal' value can be caused by finite coherence. In the context of interference 
correction two parameters are important at zero magnetic field: the thermal length 
It and the phase-breaking length l v . The thermal length, sometimes also termed the 
coherence length, is given by It = (hD/hsT) 1 / 2 . With D = 17cm 2 /s, we obtain 
It ~ 240 nm. At the lowest temperature, the devices crosses therefore over into the 
fully coherent regime. Hence, the reduced value of SG must primarily originate from a 
finite phase-breaking length l v < L. 

To compare further with theory, a value of 5G = 0.69^W/Le 2 /h was predicted for 
a graphene ribbon with normal metal contacts [TU] . For the geometry of our device this 
would translate into SG ~ O.Qe 2 /h. Our estimated normal state SG value is two times 
smaller. This can be used to estimate an effective coherence length l v by using standard 
averaging along the length of the ribbon [26| , yielding G ~ 0.69^J\V/ L^L^/L) 3 ^ 2 e 2 /h. 
With L being 225 nm, we estimate a phase-breaking length of L v « 80 nm at 230 mK 
and consequently a phase-breaking time t v « 4ps. 

The bias dependence of SG has not been studied systematically in graphene-super- 
conductor devices. However some data are available in the literature. In the recent 
experiments on single graphene layer contacted with Pt/Ta superconducting leads the 
conductance variance at 60 mK was found to be SGgs = 2Ae 2 /h [8]. Though this value 
is substantially larger than ours, a comparison has to take the larger width W into 
account (W = 2.7 /xm and L = 330 nm). It turns out that also in this experiment the 
measured value is lower than the full coherence value, which we estimate to 4.3e 2 /h 
using the expression 2.07 x 0.69^JW/Le 2 /h. 

In nanostructures made of InAs nanowires contacted with Ti/Al leads [27J ESI EH] 
the saturation value of the conductance variance in the superconducting state was found 
to be ~ 0.8e 2 /h at 22 mK J27J [28] and ~ 0A7e 2 /h in the normal state yielding an 
enhancement of ~ 1.6, whereas in |29j at 300 mK SG ~ 0.7e 2 /h in the superconducting 
state was found with an enhancement factor of ~ 1.5 compared to the normal state. 
Both the normal state values and the enhancement factors for the superconducting state 
of these results are in good agreement with our observations in a Ti/Al contacted few 
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layer graphene nanoribbon. 

In recent experiments with Niobium contacted InAs nanowires a surprisingly large 
enhancement factor of ~ 30 was found [29]. The reason for such a strong enhancement 
is that the coherence length of the nanowire and the dephasing length l v is much 
larger than the distance between the contacts, so that multiple Andreev charge transfer 
which contributes to fluctuations can occur at the interface. The increased conductance 
variance signals the transition into the superconducting state of the whole device. 

4. Summary and discussions 

In conclusion, our measurements of a few-layer graphene nanoribbon contacted with 
Ti/Al leads show pronounced UCF-type conductance fluctuations. We observe 
a decrease of the conductance variance with applied source-drain voltage with a 
characteristic cross over for bias voltages corresponding to the superconducting energy 
gap A. For voltages below A/e the conductance variance is enhanced by a factor ranging 
between 1.4 and 1.8 which is close to the theoretically predicted value of 2.07. The finite 
phase-breaking length L v < L at the base measurement temperature of 230 mK is the 
reason for the remaining discrepancy. 
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